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Abstract 

We consider an extended Hubbard model of interacting fermions on a lattice. 
The fermion kinetic energy corresponds to a tight binding Hamiltonian with near- 
est neighbour (— t) and next nearest neighbour (t 1 ) hopping matrix elements. In 
addition to the onsite Hubbard interaction (U) we also consider a nearest neigh- 
bour repulsion (V). We obtain the zero temperature phase diagram of our model 
within the Hartree-Fock approximation. We consider ground states having charge 
and spin density wave ordering as well as states with orbital antiferromagnetism 
or spin nematic order. The latter two states correspond to particle-hole binding 
with d x 2_ y 2 symmetry in the charge and spin channels respectively. For t' = 0, 
only the charge density wave and spin density wave states are energetically stable. 
For non-zero t', we find that orbital antiferromagnetism (or spin nematic) order is 
stable over a finite portion of the phase diagram at weak coupling. This region of 
stability is seen to grow with increasing values of t' . 
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The phase diagram of the 2-d Hubbard model has been a subject of intense study over 
the past few years, ever since the discovery of the high temperature superconductors. The 
strange properties of these materials are believed to arise because of strong electronic 
correlation and the 2-d Hubbard model or its strong coupling limit, the t-J model, 
have frequently been used as a starting point for theoretical attempts at understanding 
these phenomena. Even earlier, this model has been studied as the minimal model 
necessary for describing a host of phenomena ranging from metal-insulator transitions 
in the transition metal oxides to magnetism in itinerant electron systems. 

On a square lattice, at half-filling, a tight-binding Hamiltonian with hopping re- 
stricted to nearest neighbours, has a Fermi surface that is perfectly nested, so that at 
zero temperature the ground state of the system is a two sublattice spin density wave 
(SDW) state for arbitrarily small onsite repulsion (U). Detailed studies of the SDW 
Hartree-Fock (HF) solutions as well as a random-phase-approximation (RPA) study of 
the collective excitation spectrum have shown that this approach is effective even at large 
values of U where the magnetization and spin wave dispersion obtained from a mapping 
to the Heisenberg model are recovered. The introduction of a nearest-neighbour repul- 
sion (V) introduces a tendency towards charge ordering in the form of a charge density 
wave (CDW) and for sufficiently large values of V (V > Z7/4), the CDW state becomes 
energetically preferred to the SDW state. 

It has been known since the early work of Halperin and Rice that there exist 
forms of particle-hole ordering other than the CDW and SDW state. In this context, 
some other ground states which have been considered are orbital antiferromagnetism 
(OAF) and spin nematic (SN) [3-6]. These ground states can be thought of as arising 
from d-wave pairing in the particle-hole charge and spin channels respectively 0. The 
OAF and SN states have charge and spin currents in every elementary plaquette of the 
square lattice with neighbouring plaquettes having their currents oriented in opposite 
directions. These states are thus similar to the staggered flux phases which were studied 
earlier || as possible ground states of the t-J model, in the strong correlation limit. 

In this paper we present a Hartree-Fock study of the zero temperature phase diagram 
of the extended Hubbard model at half-filling on a square lattice especially focussing on 
the competition of the CDW and SDW states with the OAF or SN states. Our model has 
a kinetic energy corresponding to nearest-neighbour and next nearest-neighbour hopping 
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whose amplitudes are t and t' respectively and an onsite (U) and nearest neighbour (V) 
repulsive interactions. For if = we find that the CDW and SDW states are the only 
stable ground states of the system and are separated by a phase boundary at V — U/4. 
The SDW and CDW order parameters do not vanish at the phase boundary and the 
transition is discontinuous. However, for finite values of t' the OAF and SN states 
which are energetically degenerate within the HF approximation, are found to be lower 
in energy than the CDW and SDW states over a finite portion of the phase diagram, 
whose size grows with increasing t'. 

Our main findings (summarized in Fig. 2) indicate that the absence of nesting (which 
is destroyed by a finite t') favours the occurence of the OAF and SN states. This is 
significant since most real materials have non-nested Fermi surfaces. Stable OAF and 
SN states are found to occur in the weak correlation limit (U, V less than the bandwidth) 
where the HF approximation is expected to be reliable. The d-wave character of the 
particle-hole binding involved in the formation of the OAF and SN states, is interesting in 
the light of the observation of a pseudo-gap with a d-wave symmetry in photoemission 
experiments || [K| in the underdoped cuprates which has been interpreted variously 
as evidence for a paired state where superconductivity has been destroyed by phase 



fluctuations [11] or a manifestation of the spinon dispersion predicted by slave boson 



theories of the t-J model |]12[ . It is worth pointing out that d-wave pairing in the 
particle-hole channel, such as occurs in the OAF and SN phases, also leads to a similar 
one-particle excitation spectrum which is qualitatively consistent with the photoemission 
data. 

We consider the extended Hubbard model on a square lattice, whose Hamiltonian is 
given by 

h = Uj4, ff c jia -vY. K° + U J2 ^,t™u + y EE h i>a h i+Snni(T > (i) 

The hopping matrix elements are chosen to be = —t for i, j nearest-neighbours, 
tij = if for i, j next nearest neighbours and zero otherwise. Here /i is the chemical 
potential and U and V are the onsite and nearest neighbour repulsive interactions re- 
spectively. We work throughout at half-filling with one electron per site on the average. 
We consider mean-field states characterized by an anomalous non-zero ground state ex- 
pectation value of the average (c£ qCic+q,^) with Q = (ir, ir) being the nesting vector of 
the non- interacting Fermi surface for if — 0. 
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The CDW order parameter results from a real space decoupling of the form 

1 p 

(n it a) = 2 + 2 cos (Q- Ri ) ( 2a ) 
whereas the SDW state (we consider a z-polarized state) is described by 

1 777, 

(n iia ) = 2 + ° -J COS (Q- R ( 2b ) 

a = +(— )1 for up (down) spins. These order parameters describe a periodic modulation 
in the charge and spin density respectively resulting in a doubling of the unit cell, p and 
m are the amplitudes of the induced density modulation. 

The OAF state on the other hand has non-zero intersite averages of the form 



x Ci,a C i±xa) = ig cos(Q.Ri) (3a) 

and 



[4,a c i±y°) = -ig cos(Q.Rj) (36) 
The SN state (we consider a z-polarised state) has 

(4,aCi±xa) = i£<r cos(Q.Rj) (4a) 

and 

4 a C i±ya) = COs(Q.Ri) (46) 



It is easy to see from Eqns.(3) and (4) that the OAF and SN states have charge and spin 
currents (whose magnitude is proportional to g and I respectively) moving along the 
elementary square plaquettes of the lattice with neighbouring plaquettes having their 
currents in opposite directions (see Fig.l). These states have been discussed earlier 
by several workers [3-6]. However, their stability with respect to the CDW and SDW 
states has never been precisely established. Earlier calculations [3-5] have approximately 
evaluated the appropriate RPA susceptibilities in the normal state. The divergence of 
these susceptibilities gives the mean field transition and the temperature where this 
happens is the mean field transition temperature. It is well-known ]13[ that in two 



dimensions, enhanced fluctuations prevent long range order from developing at any finite 
temperature. The RPA fails to take proper account of these enhanced fluctuations and 
thus predicts a spurious finite temperature phase transition. The system continues to be 
a Fermi liquid below the mean field transition temperature. However, this temperature is 
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a crossover scale below which the Fermi liquid begins to develop short range correlations 
appropriate to the type of broken symmetry considered. The corresponding correlation 
length increases as the temperature is reduced and diverges at T = 0. A comparision of 
the mean field transition temperatures of the different states is a rough indication of their 
relative stability. However, an accurate determination of the zero temperature phase 
diagram of the model Hamiltonian requires a complete solution of the self-consistent 
equations (obtained by minimizing the ground state energy) for the chemical potential 
and the order parameter. A comparision of the corresponding ground state energies 
then fixes the stable ground state for the parameters considered. The resulting phase 
diagram is the main result of this paper. 

We first obtain the self-consistent equations for the order parameter and the chemical 
potential (which has to be determined self-consistently for t' ^ 0). On substituting the 
mean field ansatz of Eqns(2-4) in the Hamiltonian of Eqn(l) we find that the mean field 
Hamiltonian, for the various states considered, is given by 

RBZ 

H'MF = E [(& ~~ ^) C k,a C k,a + (£k+Q ~ ^) c k+Q,a C k+Q,a 
k,a 

RBZ 

+ E [(A^(k)c£ iaCk+Qi/3 + h.c] + X s (5) 

k,a,/3 

The reduced Brillouin zone (RBZ) is formed due to the modulation with wavevector 
Q, which introduces a two-sublattice structure resulting in a folding up of the original 
Brillouin zone. The tight-binding dispersion is £k = — 2£[cos k x -\-cosk y \-\-At' cosk x cosk y 
and A* ^(k) is the generalized order parameter where s = 1,2,3, & 4 correspond to 
CDW, SDW, OAF and SN states respectively. Thus, 

,U_ 
~2 



A* i/3 (k) = 6 atf) (— - AV)p for the CDW state (6a) 



^a,/?(k) — ° Z a,p(~~y) m f° r the z — polarized SDW state (66) 

A^g(k) = i5 at p{cos k x — cosky)(—2V)g for the OAF state (6c) 

and A^ /3 (k) = %a z a «(cos k x — cos k y )(— 2V)i for the z — polarized SN state (6c?) 

where a z is a Pauli matrix. In writing Eqn(5) we have ignored the Hartree contribution 
(U/2 + AV) which shifts the energy of each state by an equal amount and can therefore 
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be absorbed in the chemical potential. The constant piece X s in the respective ground 
state energies is given by 

X l = (2V-^)p 2 N (7a) 

X 2 = ^m 2 N (7b) 

X 3 = AVg 2 N (7c) 

X 4 = 4V£ 2 N (7d) 

where N is the total number of sites. In writing Eqn(7) we have once again ignored the 
Hartree contribution of —{U /4 + 2V)N which merely shifts the ground-state energies of 
each of the mean field solutions by an equal amount. 

The mean field Hamiltonian of Eqn(5) is quadratic and can be diagonalized by a 
canonical transformation. The folding of the Brillouin zone, caused by the (ir, it) mod- 
ulation, results in the formation of two bands whose dispersion is given by 



e^(k) = ( Vk -f l ) + V4 + A*(k) (8a) 

and 

e<"> (k) = fa -fx)- V^ + A^k) (86) 

Here A 4 (k) = (U/2 - 4V)p, A 2 (k) = -Um/2, A 3 (k) = -2Vg(cosk x - cosfc,,) and 
A 4 (k) = -2V£(cosk x - cosk y ), rj k = (^ k + £ k+Q )/2, e k = (^ k - £ k+Q )/2 and ^ k is the 
band dispersion in the non interacting limit. For t' = 0, nesting of the Fermi surface 
ensures that r^ k = 0. Then, the chemical potential stays pinned to it's non-interacting 
value (fx = 0) and the lower band (e k •*) is completely full whereas the upper band is 
completely empty. However when t' becomes finite, r^ k ^ and the bands begin to 
overlap. It is easy to see that CDW (Ai(k)) and SDW (A 2 (k)) states are isotropic (s- 
wave) whereas the OAF and SN states (A 3 (k) and A 4 (k))) states have order parameters 
of d x 2_ y 2 character. 

Using the diagonalized form of the Hamiltonian, we obtain the self-consistent equa- 
tions for the chemical potential fx and the order parameters, which are given by 



1 1 RBZ 

^ = ^j:H-^(k)) + Q{-e^(k))} (9a) 
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and 




i [e(- e (-)(k)) + e(- e (+)(k)) 



(96) 



Here Ti = 4V — U /2, T2 = U /2, Y% = Y ^ = V/2 and is the symmetry factor of the 
order parameters with Wy. = wl\ = 1 (for SDW & CDW) and w^\ = Wy\ = (cos A^ — cos /c y ) 2 
(for OAF & SN). 

We have solved Eqn(9) numerically to determine the self-consistent values of the 
chemical potential /i and the gap function A s (k). Self-consistent solutions for all the 
four mean field states exist for almost all values oiU/t and V/t and therefore it becomes 
necessary to compare the ground state energies in order to decide the relative stabilities 
of these states. The ground state energy of a given state can be easily found from the 
diagonalized Hamiltonian and is given by 



It is easy to see from a comparison of Eqns(9) and (10) that the OAF and SN 
states are degenerate within the Hartree-Fock approximation for our model. Further, 



are very different physically. The OAF state has a staggered orbital magnetic moment 
produced by the circulating charge currents which should be observable in magnetic 
neutron scattering experiments, whereas the spin-currents present in the SN state would 
not directly couple to the neutron magnetic moment. However, they are identical in 
their one-particle properties as far as experiments which probe their band structure are 
concerned. 

We use the self-consistent values of A s (k) and /i, obtained from the solution of 
Eqn(9), to evaluate the ground state energy E s in Eqn(10) for different parameter values 
t'/t, U/t and V/t. For t' = 0, we find that the CDW and SDW states are the only stable 
states and are separated by a phase boundary at U = 4V. In Fig. 2 we present the phase 
diagram obtained from Eqns(9) and (10) for typical values of t'/t in the (U/t, V/t) plane. 
Our results indicate that the destruction of nesting, introduced by t' helps to stabilize 
the OAF/SN states at weak coupling. The various order parameters do not vanish at the 
phase boundaries which therefore indicate discontinuous phase transitions. The phase 
boundaries for small values of the interaction parameters U/t and V/t have not been 
shown as numerical difficulties prevented an accurate determination of their positions. 



r) HBZ 

e s = x s + -j: [ e (-)(k) e(- e (-)(k)) + 4+)(k) e(- e «(k)) 



(10) 



the dispersion of the energy bands (e^ )(k) and e^(k)) is also identical. These states 
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The zero temperature instability of the Fermi system to the states considered by us 
can be roughly understood as arising from the energy gained due to the opening up of a 
gap A(k) in the energy spectrum which pushes down (in energy) the occupied states. It 
is then easy to see from Eqn. (6) that the gap functions in the CDW, SDW and OAF/SN 
states are proportional to 4V — U/2, U/2 and 2V respectively Thus the OAF/SN states 
are expected to be energetically close to the other two states near the U = AV line. At 
large values of U, the dominant consideration is to minimize the double occupancy in 
the ground state which is achieved by the SDW state with m — > 1. On the other hand 
for large values of V, the dominant tendency is for all the electrons to sit on one of the 
sublattices thus completely avoiding the energy cost associated with having electrons on 
neighbouring sites. This corresponds to the CDW state with p — > 1. Thus the OAF/SN 
states are expected to arise only at weak coupling and close to the U = AV line. This 
is borne out by our numerical results. The crucial role of t' in stabilizing the OAF/SN 
states is apparent from our phase diagram. The physical basis for this is not clear to 
us. However it suggests that the OAF/SN states are more robust with respect to the 
destruction of the nesting of the Fermi surface than the CDW or the SDW states. 

In Fig. 3 we show the dispersion ^ and ej^ for the OAF (SN) states as a function of 
k for t'/t finite with U/t and V/t chosen to ensure the stability of the phases considered. 
It is clear from Fig.3 that for t' ^ the upper and lower bands overlap. Thus the valence 
band is not completely full and the conduction band is not completely empty. The OAF 
(SN) state is therefore a poor metallic state with a pseudogap at the Fermi surface. 

Finally, we summarize the main results of this paper. We have obtained the zero 
temperature phase diagram of the extended Hubbard model on a square lattice at half- 
filling. We find that for finite values of next-nearest-neighbour hopping amplitude t', 
the OAF/SN states have lower ground state energy than the SDW and CDW state over 
a finite range of parameters in the weak coupling region. The size of this region, where 
OAF/SN states are stable, increases with increasing t'. 

We would like to thank A. Taraphder, R. Pandit and T. V. Ramakrishnan for useful 
discussions. 
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FIGURE CAPTIONS 

Fig.l. Schematic diagram of the circulating currents associated with the OAF and SN 
states. The arrows indicate the directions of the currents which are staggered and move 
oppositely on neighbouring plaquettes. 

Fig. 2. The mean-field phase diagram of the t — t' extended Hubbard model. All the 
phase boundaries indicate discontinuous transitions. The region of stability of the OAF 
(SN) state grows with increasing t'/t ratio, as seen in the figure. 

Fig. 3. Plot of the quasiparticle dispersion energy ej^ and e£~^ for the OAF and SN states 
along the symmetry directions in the reduced Brillouin zone. The solid curve indicates 
the conduction band dispersion while the dashed curve is for the valence band dispersion. 
The parameters choosen for the quasiparticle dispersion are t'/t = 0.1, U/t = 1.5 and 
V/t — 0.37 for which the OAF (SN) state is stable. Inset shows the position of the 
points T, X and E in the Brillouin zone. 
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